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Abstract 

We calculate the cross section for gqg-production in diffractive DIS with finite fermion masses and 
zero momentum transfer. The calculation is done in the leading log(l/a;p) approximation and is valid 
for the region of high diffractive masses (small /3) . We apply our cross section formula to diffractive 
charm production at HERA: in a Monte-Carlo-Simulation of diffractive D f± meson production we 
include both massive qq- and massive gqg-production, and we compare with preliminary HI results. 
After adjusting an infrared cutoff parameter of our model which affects the overall normalization 
of the cross section we find reasonable agreement with the data. In particular, the slope of the 
/3-distribution can be reproduced. 

1. In the process of diffractive deep inelastic scattering, 7* +p — > p+X, one can separate perturbative and 
non-perturbative contributions by filtering out particular diffractive final states. Examples of diffractive 
states which are perturbatively calculable are longitudinal vector particles or final states which consist of 
hard jets (and no soft remnant). In the latter case the hard scale which allows the use of pQCD is provided 
by the large transverse momenta of the jets, and the Pomeron exchange is modelled by the unintegrated 
gluon density. Another particularly interesting example is diffractive charm production, since the charm 
quark mass justifies pQCD, even for not so large transverse momenta of the outgoing quarks and gluons. 
Calculations for the diffractive production of massless open qq states and of massless qqg states have been 
reported in 0, g § and in @, resp., and a comparison of diffractive two-jet and three-jet events observed 
at HERA with these calculations has been presented in ||. Final states with finite quark masses have 
been calculated, so far, only for qq production Q which is expected to be the dominant final state in 
the region of small diffractive masses (large (5). However, as there are recent results from measurements 
of HI and ZEUS at HERA on diffractive production of (charmed) D*^ mesons, which extend into the 
small- /3-region, gluon radiation can certainly not be neglected, and a comparison of the perturbative two 
gluon model with HERA data has not been possible yet. In this letter we report on a calculation of 
massive qqg-production in DIS diffraction, and we present a comparison of our cross section formula with 
preliminary HI data. 

2. We will follow the study of massless qqg-production presented in Q . In particular, we again work in the 
leading-log M 2 approximation, which limits the applicability of our results to the small /3-region. Fig.[j] 
shows the notations of the process. As in [Q we restrict ourselves to zero momentum transfer, t = r 2 = 0. 
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Figure 1: Kinematics of dim-active (795 production 



As usually, Q 2 denotes the virtuality of the photon, \JW 2 the energy of the photon proton system, M 
the mass of the diffractive system, and x = Q 2 / (Q 2 + M 2 ), y = 2pq/2pl are the Bjorken scaling variables 
x and y (I is the momentum of the incoming electron). The variable (3 is defined as (3 — Q 2 /(Q 2 + M 2 ), 
and it is convenient to introduce the momentum fraction of the Pomeron xp = (Q 2 + M 2 )/(Q 2 + W 2 ). 
The limit we are interested in is defined as: 



Q 2 < M < W 2 . 

We use Sudakov variables fcj = onq' + fop + hi t (with q' = q + xp, k 2 t 
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-k 2 ), and we express the phase 

space in terms of y, Q 2 , M 2 , to 2 , t, k 2 , k 2 , with m 2 = m 2 q + k 2 , (m qq denotes the invariant mass of the 
qq-subsystem) . We obtain the following result: 
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Here 



and 



M u = J ^T{x v ,\ 2 )Tu 



(10) 



D(k) = ai(l-ai)Q 2 + k 2 +m 2 . (11) 

The function J 7 denotes the unintegrated (forward) gluon density which is connected with the usual gluon 
density g(x, Q 2 ) through: 

rQ 2 

dl 2 ^,! 2 ) = xg{x,Q 2 ). 



(12) 



The parameter a,\ is determined by the on-shell conditions for the final state particles: 
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1+ k|(k 1 + k 2) 3 ± ^ 



and it varies between and 1. The values of the momenta ki, k 2 and of m 2 decide which sign in eq. (|l2|) 
holds. 

The quark mass m q enters the calculations in two places. First, the phase space of the diffractive system 
(and so the parameter a.\ and the function S) depend upon the quark mass via the on-shell conditions for 
the outgoing particles. Secondly, the propagators of the internal fermion lines are modified by a nonzero 
quark mass which leads to changes in the matrixelements. Apart from the function D(k), cq.([ll|), which 
enters all the four 7*p-cross sections, an additional term containing the quark mass emerges in da^ j,,, 
eq.©. 

3. Using our analytic formulae, we have performed a numerical Monte Carlo calculation and compared 
with preliminary HI data on D* ± production in diffractive DIS (to_d* ~ 2.01GeV). There are (unpub- 
lished) HERA measurements on this process both from HI |7j and from ZEUS ||. Compared to other 
charmed mesons the D*^ mesons are easy to reconstruct. This makes them attractive objects for testing 
diffractive charm production. Both experiments made use of the decay channel 

D*+ - D\+ ow -> (K-n+)n+ ow (and c.c), 
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Figure 2: da/xp (sum of cc and ccg production) at different values of the cut on the gluon transverse 
momentum, k=k 2 cut and of the lower 1 - integration limit, l=l^ in (both in GeV 2 ). 



which has a branching ratio of 2.63% The kinematical regions where ZEUS and HI have taken their 
data are slightly different. We will focus on the comparison with the preliminary HI data |Q which have 
been collected throughout the years 1995-1997. The amount of data is still quite poor due to the very 
rare occurance of the D* in the considered process, and higher statistics will come from new data. 

We have implemented the ep cross section for diffractive massive cc production both from || and from 
our expression eq.(|^-[|) for the massive qqg production into the Monte-Carlo-Program RAPGAP 2.08 
JToj] . This program includes hadronization based on the Lund model. We have added the simulated cross 
sections of cc and the ccg production; the cross sections for D* + and D*~ production have also been 
added. a s was fixed at a value of 0.25. We have used the GRV NLO gluon density |llj to obtain the 
unintegrated gluon density. The boundary conditions imposed on the kinematic variables are the same 
as at HI Q. The electron and proton momenta (in the HERA system) are 27.6 GeV and 820 GeV, 
respectively, and we demand: 

0.05 < y < 0.7 
2 < Q 2 < 100 GeV 2 
x P < 0.04. 

The momentum transfer of the proton has to be limited, because in the experimental analysis also events 
with diffractive excitations of the proton have been included in the diffractive cross section: 

\t\ < 1 GeV 2 . 

Finally, there are two restriction (of more technical origin) on the transverse momenta of the D* mesons 
and on their pseudorapidity values rj = —In tan(#/2) (9 denotes the angle between the momentum of the 
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particle and the direction of the incoming proton) : 

in the ep-CMS: M^**)! < 1-5 and priP**) > 2 GeV. 

In addition to these restrictions which come from the experimental side there is one further cut on the 
ccg final state due to non-perturbative elements in the calculation. Since the transverse momentum of 
the final state gluon defines the scale for the strong coupling constant we have to impose a lower cutoff 
\f-2cut on the gluon transverse momentum. If it is choosen to be too large, we loose on the cross section; at 
too small values the reliability of perturbation theory becomes weak. In our calculations we have choosen 
k| cui = 2 GeV 2 . Finally, the 1 integration of the gluon loop in cq.dHj) is not defined in the infrared 
region. We therefore choose a nonzero lower integration limit l 2 nin at = 0.5 GeV 2 . In the numerical 
treatment of the cc production cross section neither k| cut nor l 2 nin are needed. 
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Figure 3: The cross section of diffractive D* production in the 2-gluon-model as a function of xp. The 
data points are preliminary HI data H. 



In order to test the sensitivity of our results to the choice of the two cutoff parameters we have computed 
da/dxp for different values. Results are shown in fig.^: the main change is in the overall normalization. 
As expected, lowering any of the two cutoffs leads to an increase of the cross section. In our subsequent 
analysis we have choosen to adjust our xp distribution (cc and ccg) of the D* mesons to the data in the 
lower xp-bin. Since there is a certain arbitrariness in this choice, in our comparison with preliminary HI 
data we will concentrate on the shape of the D* distributions rather than its normalization. 
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Turning now to the results of our numerical analysis we start with the xp distribution. The contributions 
due to ccg and cc and the sum of both are shown in fig.|3|. (as we have said before, the value of the 
sum of both cross sections (cc and ccg) has been adjusted to the data in the lower xp-bin). In the 
upper bin the computed cross section is by a factor of about 10 smaller than the data point. In this xp- 
region it is expected that, apart from Pomeron exchange (which, in our model, is the 2-gluon exchange) 
secondary exchanges have to be included: in a perturbative description such an exchange corresponds to 
qq-exchange. Since such a contribution has not yet been included, it is not surprising that the two-gluon 
model undershoots the data. 
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Figure 4: The cross section of diffractive D* production in the 2-gluon-model as a function of /3. Data 
points are preliminary HI data j7j. 



The D* cross section as a function of (3 is shown in fig.[4|. The measured cross section is rising with 
decreasing (3. This property is not reproduced by the cc channel alone. Only when both cc and ccg 
production are included the (3 distribution of the two gluon model increases with decreasing /?. This fact 
is nearly independent of the choice of the cutoffs \^ nin and k^ ut which mainly affect the overall strength 
of the cross section. Therefore the ccg contribution really improves the description of the /3 distribution. 
The calculated sum of the cross sections overshoots the data in the upper /3-bin and still remains below 
data in the lower /3-bin. A reason for the disagreement in the lower /3-bin could be the radiation of more 
than one gluon which might be present in the data but is not included in our model. Since in the small 
/3-region xp gets larger, the disagreement could also be due to our neglect of secondary exchanges. 
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Two further distributions are shown in ng.||. In contrast to the cc contributions alone the full cross 
sections in both bins lie close to the data points (except for the low-?/ point in fig.pp even within the error 
bars of the data points). The ratio of the data values in the lower and upper Qr bins is about 15. For 
cc alone the model gives the value 8.3; when the ccg contribution is included it goes up to 10.7. Also in 
the y distribution (fig.||b) the agreement of the sum of both cross sections with data is better than for cc 
alone. 
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Figure 5: The cross section of diffractive D* production in the 2-gluon- model as a function of Q 2 and as 
a function of y. Data points are the HI preliminary data Q|. 



4. In this letter we have analysed, within the perturbative two-gluon model, DIS diffractive charm 
production (production of D*^ mesons). Compared to an earlier attempt where only cc production 
had been included in the theoretical analysis the present analysis contains, as the new ingredient, also 
(massive) ccg production and leads to a considerable improvement in the agreement with experimental 
data. 

Despite this encouraging success, several improvements in the theoretical part should be made. First, 
the cross section formula for ccg production has been calculated in the leading log-M 2 approximation; 
an improvement which extends the applicability down to small-M 2 values would be very desirable. For 
consistency reasons, one then will need a NLO-calculation of cc production. Next, our comparison 
with data indicates the need of cegg final states: such an extension (at least in the leading log-Af 2 
approximation) should be fairly straightforward. Finallly, the region of xp > 0.02 seems to require 
secondary exchanges which, in the framework of perturbative QCD, should be modelled by qq exchange. 
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A successful test of the two-gluon model in DIS Diffraction, apart form providing a description of charm 
or jet production at HERA, is also of general theoretical interest: the cross section formula for diffractive 
qq + ng production contains the perturbative triple Pomeron vertex which is expected to play a vital role 
in the unitarization of the BFKL approximation. It has been calculated both analytically and numerically, 
and these calculatutions can be tested experimentally in DIS diffraction dissociation. 
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